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Abstract 

Two- and three-pion correlation functions are investigated for a source that 
is not fully chaotic. Various models are examined to describe the source. 
The chaoticity and weight factor are evaluated in each model as measures of 
the strength of correlations and compared to experimental results. A new 
measure of three-pion correlation is also suggested. 
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I. INTRODUCTION 



Two-pion correlations obtained in relativistic heavy-ion collisions have been used to 
extract the size and shape of the pion-emitting source, based on the Hanbury-Brown Twiss 
(HBT) effect. Correlations are not, of course, limited to those of two pions, but can also be 
of multi-pions. Though multi-pion correlations are unavoidably complicated, they contain 
new information not available from the two-pion correlations. For example, when final-state 
interactions of the emitted pions are neglected, the two-pion correlations for a chaotic source 
depend only on the magnitude of the Fourier transform of the pion-source function, but the 
multi-pion correlations depend also on its phase |]T],0]. The magnitude is an even function of 
the relative momentum of the emitted pions, while the phase is an odd function. One thus 
hopes to extract new information about the source from the multi-pion correlations. 

In a previous article 0, we made a detailed investigation of three-pion correlations, 
the simplest multi-pion correlations after the two-pion correlations, over a wide range of 
kinematics, to investigate the feasibility of extracting new information about the source 
through the phase in the case of a chaotic source. We find that extracting new information 
would be rather difficult in practice because the multiplicative factor of the phase function 
becomes small in the region where the phase varies appreciably. 

In experiments, the two-pion correlations at the zero relative momentum is observed to 
be less than two. After the final-state interactions are removed, it must be two in the case 
of a chaotic source. The measurement suggests that the source is not completely chaotic. 
Recently, a measurement of the three-pion correlations has been reported, showing that the 
strength of the three-pion correlations are also less than expected in the case of a chaotic 
source. 

In this paper we investigate the two- and three-pion correlations for a source not com- 
pletely chaotic. We introduce models of the source that are various mixtures of coherent and 
chaotic sources. The models also include one with a novel structure, a mixture of multiple 
coherent sources and a chaotic source. These models describe different dynamics generating 
pion emissions, though we do not pursue the identification of dynamics corresponding to 
each model. 

A mixture of a coherent source and a chaotic source is a model that has often appeared in 
the literature under the name "partially coherent source." We find that this popular 

model poorly reproduces the recent data of the two- and three-pion correlations ||. A 
model consisting of multi-coherent sources and one chaotic source appears to yield a good 
agreement. 

It is both theoretical and experimental practice to calculate the so-called weight factor 
from the three-pion correlations at the vanishing relative momenta. The weight factor is 
usually considered to represent the strength of the genuine three-pion correlations. When we 
examine more complicated models than the popular partially coherent model, we find that 
the weight factor no longer describes the genuine strength. The expression that yields the 
genuine strength depends on the structure of the source, and there is no universal expression 
as such. Nevertheless, we propose a modified expression for the weight factor that has a 
wider validity. 

In Sec. U we define pion spectra, correlation functions, and measures of two- and 
three-pion correlation functions, chaoticity, and weight factor, respectively. In Sec. Ill, 
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various nonchaotic models of the pion-emitting source are introduced, and chaoticities and 
weight factors are obtained for them. Section IV presents discussions and a summary: a 
new expression of the weight factor is introduced and discussed, and chaoticities and weight 
factors are compared with the recent experimental data. Appendix [A] gives a derivation of 
correlations for multi-coherent sources. In Appendix [B|, the new weight factor is derived. 



II. CORRELATION FUNCTIONS, CHAOTICITY, AND WEIGHT FACTOR 

For the sake of clarity in the sections to follow, we define correlation functions and their 
measures, chaoticity, and weight factor. Note that our definition is standard. 

We first write the basic quantities, pion spectra, in second-quantization form, as follows: 

Wi(p) = (ala p ), (1) 
W 2 (p 1 ,p 2 ) = (a P i a P2 a pi a P2)> ( 2 ) 
W 3 (p 1 ,p 2j p 3 ) = (a P i a P2 a P3 a Pi a P2 a P3)- ( 3 ) 



Though these expressions are simple and reasonable, the explicit definition of (...) is a 
complicated issue. (...) represents a quantum statistical average and is formally written as 
in terms of the quantum state, or Tr{/S...} in terms of the the density matrix, 
p. In this work, we introduce models to represent various underlying dynamics of the pion 
emission. Note that the momenta above and hereafter are on-shell, e.g., p° = y/p' 2 + m' 2 for 
p. 

In terms of the spectra, we define two- and three-pion correlation functions in the usual 
way P|, 

n ( \ (n) 2 W 2 ( Pl ,p 2 ) 

C 2 {p u p 2 ) = — — , 4 

(n(n - 1)) VKi(pi)Wi(p 2 ) 

3\Pi,P2,P3) ( n ( n - i)( n - 2 )) W 1 ( Pl )W 1 { P2 )W 1 ( P3 y 1 ' 

Here, the normalizations are introduced in order to take account of multiplicity fluctuation, 
with the following definitions: 



(n> = j rfWi(p), (6) 
(n(n - 1)) = J d 3 Pl d 3 p 2 W 2 { Pl ,p 2 ), (7) 
(n(n - l)(n - 2)) = f d^d^d^W^p^ps). (8) 



In the usual models, the above normalizations yield C 2 and C 3 — > 1 as the relative momenta 
approach infinity. In more complicated models, such as those that possess the particle- 
number fluctuation per mode, however, the normalizations of Eqs. @ and (|5[) yield an 
asymptotic value different from unity. We will discuss this point fully in Sec. [1 1 j~D| . 

When the HBT effect does not appear, correlation functions are independent of relative 
momenta. When the HBT effect occurs, however, the correlation functions are no longer 



3 



constant, but are dependent on relative momenta. As is well known, the size and shape 
of a pion-emitting source are extracted from the relative-momentum dependence of the 
correlations through the HBT effect. The methods of extracting the size and shape of the 
source has been discussed numerous times, and we will not go into the issues here. 

We define the chaoticity, X(p), and the weight factor, cj(p), which are usually considered 
to represent measures of the strength of the two- and three-pion correlations, respectively 



\(p) = C 2 (p,p)-l, (9) 

/ \ C 3 (p u p 2 ,p 2 ) - 1 - (C 2 (pi,p 2 ) - 1) - (C 2 (p 2 ,p 3 ) - 1) - (C 2 (p 3 ,pi) - 1) 
u(p) = 



2^2(^2) - 1)(C 2 (P2,P3) " l)(C 2 (p 3 ,Pl) " I] 



Pl=P2=P3=P 

(10) 



Note that we define the weight factor explicitly at pi = p 2 = p 3 , but that it has been 
extracted from measurements over small relative momenta |J. The two methods show 
practically no difference since the weight factor is expected to be a slowly varying function 
of the relative momenta . 

The weight factor, however, is not always the measure of the strength of the genuine 
three-pion correlations. Generally, the three-pion correlation functions are related to the 
two-pion correlation functions as 

Cs{pi,P2,P3) = 1 + a(C 2 (p!,p 2 ) + C 2 {p 2 ,p 3 ) + C 2 {p 3 ,pi) - 3) 
+ (genuine three-pion correlation) 

+ (other two-pion correlation), (11) 

where the coefficient, a, is not always unity even at p± = p 2 — p 3 . Consequently, Eq. (]To|) 
does not always remove the linear dependence of C 2 s from C 3 . In the cases of a chaotic 
source and a partially coherent source, a is unity, but for sources of more complicated 
structure, it is not. Furthermore, C 3 generally depends on the two-pion correlations not 
only as the linear C 2 s but also in more complicated ways as a function of two momenta. 
We will discuss these points fully in the following section. 

The relation between C 3 and C 2 depends on dynamics and is model-dependent in practice. 
The extraction of the strength of the genuine three-pion correlations thus requires explicit 
knowledge of the dynamics. Accordingly, in this paper we adopt the conventional approach 
of using the weight factor defined as Eq. (0), for numerical results as a convenient means 
for making comparisons among various models. 

III. DIFFERENT MODELS OF THE SOURCE 

A. Coherent Source and Chaotic Source 

The spectra are often written by c-number source current [f3]||, yielding a chaotic source 
and a coherent source as the two extreme cases. The chaotic source shows the HBT effect 
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with both the chaoticity and the weight factor as unity. The correlation functions for a 
coherent source are independent of the relative momenta, and the chaoticity for it is zero. 

The difference between the out-state and in-state annihilation operators of an emitted 
pion defines the source current as 

flout (p) = flin(p) +i [ d^x 1 =J(x)e- ip ' x 
J ■ /(2vr) 3 • 2p° 



= a in (p) + i 1 J(p), (12) 
^(2vr)3 • 2j9° 

where J(x) and J(p) are the source current and its Fourier transform, respectively. The 
average that appears in Eqs. (J]])-© can be written as 

(•••) = / vr(p)vj( P )v[j*(p), JipWinl ■ ■ ■ \0in), (13) 

where |0j n ) is the in-state vacuum and V[J*(p), J(p)] is the distribution functional of J(p), 
which has a statistical fluctuation. The pion spectra are obtained in this case from Eqs. 



^ ) = (27F^ (|J(p)|2)j ' (14) 

W 2 ( Pl ,p 2 ) = {2n y] p o. 2p o (\ J (Pi)\ 2 \AP2)\ 2 )j, (15) 
Ws(PnP2,P3) = {2n)9 2p o. 2p o . 2p o (\APi)W(P2)W(P3)\ 2 )j, (16) 

where (• • -)j = / VJ*(p)VJ(p)V[J*(p), J(p)] • ■ ■■ Hereafter, we do not explicitly show the 
subscript J of the angle bracket. If the phase of the source current is random, we call 
the source chaotic. In this case, V[J*(p), J(p)] have a Gaussian form and higher-order 
moments of J(p) are represented by the second-order moment such as 

(J*(pi)J*(p2)J(qi)J(q2)) = 

(J*(ft)J(ft))TOJ(ft)) + (r( Pl )J(q 2 ))(r(P2)J(qi))- (17) 

For a chaotic source, the pion spectra are written as 

Wi(pi) = F n , (18) 
W 2 (pi,p 2 ) = F n F 22 + F 12 2 , (19) 
W 3 ( Pl ,p 2 ,p 3 ) = F n F 22 F 33 + F aaF bc 2 + 2F 12 F 23 F 31 cos(<l> 12 + $ 23 + $ 3 i), (20) 

(a,b,c) 

where the amplitude, F a b, and the phase, $ a {, are defined as |l[ 

F o6 exp(^ a6 ) = / A J*(Pa)J(pt)). (21) 

(2tt)3 J2pl ■ 2pl 
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S(a,6,c) implies a sum over (a, b, c) = (1, 2, 3), (2, 3, 1), (3, 1, 2). The correlation functions are 
then 

C 2 ( Pl ,p 2 ) = l + -p^, (22) 

^11-^22 

C 3 ( Pl , P2 , P3 ) = 1 + E + 2 l 12 l 23 l 31 COS ( $ 12 + $ 23 + $81), (23) 

(a,b) ^aa-Fbb -^11 -^22 -^33 

where Z)( a j) is over (a, 6) = (1, 2), (2, 3), (3, 1). F ab and $ ab contain the information about 
the size and shape of the pion-emitting source, due to the HBT effect. The chaoticity and 
weight factor for a chaotic source always achieve unity. 

When the source current has no randomness, the source is called coherent. In this case, 
the n-pion spectra are expressed as 

n 

W n (p 1 ,..., Pn ) = l[W l (p l ). (24) 
i=i 

The HBT effect does not appear in this case. The correlation functions achieve unity, and 
the chaoticity vanishes. 



B. Partially Coherent Source 

In nuclear collisions, one may expect to involve the mixture of a coherent source and a 
chaotic source, as suggested by the observation that the chaoticities are often smaller than 
unity. Here, we sketch the case of the partially coherent source, mostly following Ref. ||T|. 

When the pion-emitting source is a mixture of a chaotic source and a coherent source, 
the source current is written as a sum of both currents, J(x) = J c ha(%) + Jcoh(x) ULi]- The 
sources are not correlated with each other, (J* ha (p)J co h(q))=0. The one-pion spectrum and 
correlation functions for this source are obtained as 

Wi(pi) = /ii + Fn, (25) 

n , \ -i . F 12 2 + 2/iaFia CQS($1 2 - 012) 

C2iPl > P2) = 1 + (/n + F 11 )(/22 + F 2 2) ' (26) 

,1 , \ 1 , V- F <* + 2fabF ab COs(® ab - (j) ab ) 

C 3 (PI,P 2 ,P3) = 1+2^ 



(a,b) 



( faa F aa ) ( f bb + F t 



bb) 



where 



+2 (/„ + F„)(fa!^)(/ 3 3 + F 3 3) { f '^ A C0S< * 12 + * 23 + * 3l) 
+ E lfabF bc F ca cos((P ab + $ 6c + $«,)] |, (27) 

(o,fe,c) 



f ab exp(i(j) ab ) = Jcoh(Pa)Jcoh(Pb), (28) 

(2ic)^J2jfi a ■ 2pl 
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and F a b and $ a fc are the same as those for a chaotic source, Eq. (|2T|). The chaoticity and 
weight factor are evaluated as 

A(p) = e(p)(2-e(p)), (29) 



where e(p) is the fractional parameter of the coherent source: 

6(Pl) = T^V- (31) 
ill + -Til 

The source becomes chaotic for e(p) = 1, while it becomes coherent for e(p) = 0. Figure [I] 
shows the chaoticity and weight factor as functions of e(p), and in Fig. |^ the weight factor 
is shown as a function of chaoticity, varying as e(p). The weight factor can vary from to 
1. In this case, the term — 3A, in Eq. (|TUp , completely removes the two-pion correlations. 



C. Multi-Coherent Sources 

We consider the mixture of a small number of coherent sources. This model differs from 
the case of multiple coherent sources that was previously examined in and, to the best of 
our knowledge, this is the first time that this model has been explicitly discussed. We assume 
here that the pion-emitting source is made of N coherent sources that are not coherent with 
each other and appear to be obeying the Poisson distribution. Note that the model considered 
in [[3] is a mixture of multiple coherent sources that are randomly distributed, and that the 
model was introduced as a description of a chaotic source when the number of the coherent 
sources becomes large. In a sense our model here is the opposite limit of the small number 
of coherent sources in the model considered in 0. In the following, we present merely the 
final expressions, leaving their derivations for Appendix A. 

The source current can be written as 

N 

J(x) = Y,j(x-X n )e- i0 ", (32) 

n=l 

where j(x — X n ) is the n-th coherent source current, located at X n with the random phase 
6 n . Each coherent source is assumed to be expressed by the same j(x) but to be located at 
a different position. There are N coherent sources, obeying the Poisson distribution. Note 
that the distribution must be renormalized in order to exclude the no-source event that is 
not observed. 

The one-pion spectrum and normalized correlation functions are obtained as 

W( P i) = a!—^, (33) 
1 — e a 

C 2 ( Pl ,p 2 ) = l + ^—\p l2 \ 2 , (34) 
a + 1 

oc(cn + 2) o? 
Cs(p hP2 , P3 ) = 1 + a2 + 3a + 1 £ IP^I 2 + a 2 + 3a + 1 ■ 2Re (^3P3l) , (35) 

(a,b) 
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where pij is the Fourier transform of the spatial distribution of coherent sources, and a is 
a parameter of the Poisson distribution. The mean number of coherent sources is a/(l — 
exp(— a)). In this case, the subtraction — 3A, in Eq. (Jl0|), does not remove the two-pion 
correlations completely, because the coefficient of \pi 2 \ 2 in Eq. fl34|) is different from that in 
Eq. (|35|). The chaoticity and weight factor are derived as 

A = ~~7i (36) 
a + 1 



12a 2 + 2a + 3 /a + 1 
W= 2a 2 + 3a + l V~' ( } 

In Fig. |2], the chaoticity and weight factor are illustrated as functions of the mean number 
of sources, and Fig. |^ shows the weight factor as a function of the chaoticity. The chaoticity 
varies from to 1, because this source becomes coherent at a = and chaotic at a — > oo. 
The weight factor diverges at a — > 0. The reason for this is the failure of the subtraction 
—3 A, in Eq. (|TUD, and the divergence is not caused by the genuine three-pion correlation. 

D. Multi-Coherent Sources and One Chaotic Source 

We now consider the source in which the multi-coherent sources of the previous subsection 
are mixed with a chaotic source. The multi-coherent sources are not coherent with each 
other, as before. The source current is written as 

N 

J(x) = "£j{x - X n )e- ie " + J cha (x), (38) 

n=l 

where j(x) and J c / ia (x) are a coherent source current and the chaotic source current, respec- 
tively. The n-th coherent source is located at X n and is distributed with p(X n ). The phase, 
9 n , is randomly distributed between and 2ir. There are N coherent sources, obeying the 
Poisson distribution, Eq. ( |A9| ) . 
Pion spectra come out to be 

W 1 (p 1 ) = ah u + F n , (39) 
W 2 ( Pl ,p 2 ) = W 1 (p 1 )W 1 (p 2 ) + ah n h 22 + \ah 12 e^ + F 12 e^ , (40) 
^3(^1,^2,^3) = W 1 {p 1 )W 1 {p2)W 1 {p 3 ) + (3a + l)ah u h 22 h 33 + a ^ h aa h bb F cc 

(a,6,c) 

(a,6,c) 

+2 h aa (a 2 \h bc \ 2 + ah bc F bc cos($ bc - ip, 

(a,b,c) 

+2Re {{ah l2 e^ + F l2 e^){ah 23 e 1 ^ + F 23 e i ^ 3 )(ah 31 ^ 31 + F 31 e i<I>;i1 )} , 

(41) 

where a is the parameter in the Poisson distribution, and 



-be 
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3*{pi)j(Pj) 



(2 7 r) 3 v /2rf • 2p° 2 



Piy 



(42) 



For \pi - p 2 \ 



oo, we have 



e(pO)(l-e(P2)), 



Wi(pi)Wi(p 2 ) 

where the fractional parameter of the chaotic source, e, is defined as 



e(Pi 



a/in + Fu 



(43) 



(44) 



The normalization factor of C2(pi,P2), shown in Eq. @, does not yield the proper asymp- 
totic value of unity as the relative momentum approaches infinity. Generally e depends on 
the momentum, and the asymptotic value of C 2 thus depends on the two momenta sepa- 
rately. If we assume that e is independent of the momentum, we find that C 2 is normalized 
properly by the use of Eqs. (f|) and (|5]). In this case, the correlation functions are given by 

ah 12 e^ 12 + F 12 e^ 12 2 



C 2 (P1,P2) 
C 3 (P1,P2,P3) 



1 + 



1 + 



wmwm (i + J(i 



i + ^(i-e) 3 + !(i-6) 2 6 



(45) 



E 

(o,6) 



ah ab e^ + F ab e 



W 1 (p a )W 1 { Pb ) 



2Re {(ah 



#12 + F l2 e l ^ 12 ){ah 23 e^ 23 + F 23 e^ 23 )(ah 31 e^ 31 + F 31 e i<s>31 



)} 



+2£ 

(a,b) 



a(l — e)h 



W l (p 1 )W 1 (p 2 )W 1 (p 3 ) 
(1 - e)h ab F ab cos($ afe 



ab 



W 1 (p a )W 1 (p b 



W 1 (p a )W 1 (p b ) 



(46) 



Comparison of the preceding C 2 and the first term in the square bracket in C 3 shows that the 
two-pion correlations are not to be removed completely, as a ^ 1 in Eq. ( |TTD in this case. The 
conventional weight factor thus no longer represents the strength of the genuine three-pion 
correlations. Furthermore, the last sum in the square bracket of the above C 3 represents 
the effects of the two-pion correlation (because it depends on a pair of the momenta), 
corresponding to the "other two-pion correlation" in Eq. QTTj) . 
The chaoticity and weight factor are obtained as 



A 



OJ 



a 



a + (1-e) 



2a 2 + 2a(l-e) 2 + 3(l-e) 3 (l-2e) 



la 



(1- 



2(a 2 + 3a(l-e) 



a 



(47) 
(48) 



Figures |3] and f| show A and u>, respectively, as functions of a at e = 0.1 — 0.9. In Fig. ||, 
the weight factor is shown as a function of the chaoticity, varying e for various a's. The 
divergence of u at a = 0, except for e = 0.5, is caused by the incomplete cancellation of the 
C 2 's. 
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IV. DISCUSSIONS AND SUMMARY 



We have examined various models of the source that are not completely chaotic. The 
chaoticity comes out to be between and 1 in all models, but the weight factor takes a 
wide range of the value. The value of the weight factor even diverges in some cases, as a 
consequence of incomplete removal of the two-pion correlations from C3. For a partially 
coherent source that is commonly examined, the removal is complete and the anomalous 
behavior does not appear. 

There is a way to avoid such incomplete removal, at least in all the models that we have 
examined here. It is done through a new subtraction, 

R , \ nt \ 1 W (2(n(n - 1)) - (n) 2 ) (n(n - 1)) ^ 

R3\Pl,P2,P3) = C 3 (p 1 ,p 2 ,p 3 ) - 1 — 2^ {C 2 {Pa,Pb) ~ 1) , 

[n{n L){n Z)) {n) 

(49) 

where (n) and the similar are defined in Eqs. (0) - ©• This subtraction works correctly 
for the multi-coherent sources and also for the mixture of multi-coherent sources and one 
chaotic source (but at the zero relative momenta). We find that the terms corresponding to 
the "other two-pion correlations" also vanish at the zero relative momenta (see Appendix 
B). Using Eq. ([49]) , we define a new weight factor as 

_ 3 

. (n(7i-l)(n-2))l , ,/(fi(n-l))„ f „, ,\ 5 

= (n)3 ~2 R ^P^P^ 1 (w> 2 - E{^^)-1}1 ■ (50) 

For the partially coherent source, the new weight factor becomes the conventional weight 
factor of Eq. 

Though the new weight factor properly represents the strength of the genuine three-pion 
correlations for all models examined in this work, the new weight factor must not have this 
property at all times. The relation between CVs and C3 is generally complicated, depending 
on dynamics in each case. We thus do not expect that the new weight factor has a universal 
application. 

Figures [| and ^illustrate the weight factors as functions of the chaoticities for the different 
models, together with the experimental data from the CERN NA44 Collaboration ||. In 
this experiment, the chaoticity and weight factor are measured as 0.4 — 0.5 and 0.20 ± 0.19, 
respectively. Figure § shows that the partially coherent source and the mult i- coherent sources 
disagree with the experiment. The mixture of multi-coherent sources and one chaotic source 
reproduces the data if we set a = 0.13 and e = 0.60. This corresponds to the mean number 
of the coherent sources being 0.13 and to about 60 % of the total pions emitted from the 
chaotic source. The ratio of the pion number emitted from one coherent source to that from 
one chaotic source, hn/F n , is then about 5. This ratio may be unrealistically large, but 
we note that the experimental data are "minimum-bias," suggesting that our best-fit may 
not be unrealistic since the multiplicity fluctuation can be large. Further data are needed 
to confirm that this is indeed the case. 

In summary, we investigate the two- and three-pion correlations for various models of 
a source that is not completely chaotic. The chaoticity and weight factor are evaluated as 
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measures of two- and three-pion correlations. The chaoticity always varies between and 
1, but the weight factor takes the value of a wide range and sometimes even diverges. The 
conventional weight factor includes the effects of the two-pion correlations in some models, 
yielding the anomalous behavior. We propose, in all model considered here, a new weight 
factor that has no such difficulty, but not expected to be valid universally. We find that the 
model of multi-coherent sources and one chaotic source could reproduce the chaoticity and 
the weight factor observed in the recent experiment. 
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APPENDIX A: MULTI- COHERENT SOURCE MODEL 



We assume that iV coherent sources are created during a collision and that the position 
of the n-th source is X n . The sources are uncorrected with each other. A source current is 
defined as 



N 



J(x) = J23(x-Xn)e- tdn . 



(Al) 



n=l 



where 9 n is a unique random number varying form to 2n. The average about 9 n is denoted 
by (•••)• The Fourier transform of the source current is written as 



N 



(A2) 



71=1 



where j(p) is the Fourier transform of j(x). The one-pion spectrum is 

1 



W N ( Pl ) 



(2vr)3p? 
\j(pi)\ 



o(\J(Pi)\ 2 ) 



(27r)3 p o 
ai N, 



2 N 



ipi-(X„— X m ) 



'n w m \ 



where p° is on-shell and 



ai = 



|j(Pi)| 2 
(2ti)Yi' 



(A3) 



(A4) 



The two-pion spectrum is 

W N (p 1 ,p 2 ) 



(2n)YiP°2 



a\a 2 



N 



jy2 _|_ e -i(pi-P2)-(x m -x n ) 



(A5) 



The location of each small source, X n , is assumed to obey a distribution p(X n ), which is 
normalized to be unity, or / p(X n )dX n = 1. Averaging the two-pion spectrum with the 
distribution, we obtain 



W N ( Pl ,p 2 ) = J P (X 1 )dX 1 ---p(X N )dX N W N (p 1 ,p 2 ) 



aia 2 



N 2 + N(N -l)\ Pl2 \ 



(A6) 



where p 12 = J p(x)dxe-^-^ x . Note that EijAj = £*Ai + £i#4j and 1 
AT (AT — 1). 

The three-pion spectrum is 
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W N ( Pl , P2 , P3 ) = 1 0n0 (\J(Pl)W(P2)W(P3)\ 2 }, 

(A7) 

and the averaged spectrum is evaluated as 

W N ( Pl ,p 2 ,p 3 ) = a ia2 a 3 \N 3 + N 2 (N - 1) (\p 12 \ 2 + \p 23 \ 2 + |p 3 i| 



+N(N - 1)(N - 2) • 2Re (p 12 p 23 p 3 i)] . (A8) 

Note that J2ij,k ^ijk — J2i \a + y-^-uj + -^.ji + A?'m) + J2i^j^tk ^-ijk an d that J^i^j^k 1 = 
JV(JV-l)(iV-2). 

We next assume that the coherent sources obey the Poisson distribution. The Poisson 
distribution is 

N 

P^ = ^- e - a foriV = 0~oo. (A9) 

The N = event should be excluded because in such an event no pions are emitted. The 
distribution we use is renormalized as 

Pn = ^^j foriV = l~oc. (A10) 
The expectation values change from those of the usual Poisson distribution, 

Wp = 7^ (AH) 
(N(N-1)) P = T ^, (A12) 

(N(N-l)...(N-n)) P = 1 ^. (A13) 

Pion spectra averaged with the Poisson distribution, W = J2n=i PnW n , are 

W(p 1 )=a 1I -2L_ (A14) 
W( Pl ,p 2 ) = \a(a + 1) + a 2 |p 12 | 2 l , (A15) 

I — e a i j 

W(Pi,P2,Ps) = fllfl2fl3 \ a {a 2 + 3a + 1) 
1 — e 1 

+a 2 {a + 2) (\p l2 \ 2 + \p 23 \ 2 + \p 3 i\ 2 ) + a 3 -2Re(p 12 p 23 p 3 i)] • (A16) 

The correlation functions are 

C 2 ( Pl ,p 2 ) = l + ^-\p 12 \ 2 , (A17) 
a + 1 

^ / \ ^(^ + 2) /, I o i I o | , 
<?3(Pl,P2,P3) = 1 + 2 , o rri^ 12 + IP23| + P31 

cr + 6a + 1 v 
a 2 

+ 2 ^o T 1 • 2Re (Pl2P23P3l) , (A18) 

a z + 6a + 1 
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where the correlation functions are normalized to be unity at \pi —pj\ — > oo. The chaoticity 
and weight factor are 

A = ^, (A19) 



1 2a 2 + 2a + 3 a + l 
UJ ~2a 2 + 3a + l\a' ^ ' 



APPENDIX B: DERIVATION OF EQ. (|o|) 

The two- and three-pion spectra can be written generally as 

W 2 (pi,p 2 ) = (a\ai)(ala 2 ) + (a{a 2 ) (4«i) + ((44 a i fl 2)), (Bl) 
W 3 (pi,p 2 ,p 3 ) = (a\ai}(ala 2 }(ala 3 } + ( a l a a) ( a l a c) (4 a &) + 2Re ((ala 2 )(4«3)(4 a i)) 

(a,6,c) 

+ 2 (4 a a)((44 a b a c)) + 2Re (4 a &)((44 a a«c)) + ((444 a l a 2fl3)), 

(a,6,c) (a,b,c) 

(B2) 

where aj is an annihilation operator of momentum pi. We assume that (a«) = 0, and there- 
fore we don't consider the coherent source and partially coherent source in this appendix. 
((• • •)) corresponds to a cumulant, or a connected Green function. If we have a generating 
functional, Z[z*(p), z(p)], such as 



(Q t Q t Qi , n) _ 8 2n Z[z*(p),z(p)} 

1 " ' $ Z {Pl) ■ ■ ■ Sz(p n )Sz*(p 1 ) ■ • • 5z(p n ) 

then the cumulant is obtained as 

5 2n (In Z[z*(p),z(p)\) 



2=0 



((4 • ■ ■ 4°1 ■ ' ' a n)) 



(B3) 



(B4) 



2=0 



5z(pi) ■ ■ •5z(p n )5z*(pi) ■ ■■5z(p n ) 
We introduce two assumptions. The first is 

((44 a i a 2)) = A(a\ai) (a 2 a 2 ) (B5) 
((444 a i a 2a3)) = B{a\a 1 ){a i 2 a 2 ){a\a 3 ), (B6) 

where A and B are constant. The second assumption is that the interference terms, such as 
(a\a 2 ) (4°i) > vanish by integrating over the momenta, for example, 

J d 3 p 1 d 3 p 2 (a\a 2 )(a 2 ai) = 0. (B7) 

Thus we can obtain 

(n(n-l)) = (1 + A)(n) 2 , (B8) 

(n(n-l)(n-2)) = (1 + 3A + B)(n) 3 . (B9) 
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Correlation functions at the zero relative momenta become 



C 2 (p,p) = l 
C 3 (p } p,p) = 1 + 



1 + A'' 



1 + 3A + B 



{3 + 2 + 6A} . 



(BIO) 
(Bll) 



The terms 3, 2 and QA between the braces in the above equation correspond to 
J2(a,b,c)( a i a ^( a l a c}( a l a b), 2Re ((afas) (aja 3 ) (4 a i)) and 2Re J2(aM)( a l> a ^(( a t a t a o. a c}) in Eq. 
( p2|) . To obtain the genuine three-pion correlation, we should subtract the first and third 



terms between the braces, 
discussion, 



Thus we define the new weight factor, considering the above 



as 



u — - 



1(1 

2 



3 A + B) (C 3 {p,p,p) - 1) - 3(1 + 2 A) (C 2 {p,p) 



[(1 + A) (C 2 (p )P )-l)]i 
1 <n(w - ( g"" 2)) (C a (p,P,p) - 1) - 3 • "W^-W (C a (p,p) - 1) 



^#(C 2 (p )P ) 



1) 



(B12) 



This equation corresponds to Eq. (|50|). 

The two assumptions are valid for multi-coherent sources and one chaotic source with 
setting A = a _1 (l — e) 2 and B = a~ 2 (l — e) 3 . For multi-coherent sources, the assumptions 
are not valid because of the factor (1 — exp(— a))" 1 , in Eqs. ( |A14j )- (|A16|) , due to the 
renormalized Poisson distribution, Eq. ( |A10[ ). In using the usual Poisson distribution instead 
of the renormalized one, the factor vanishes, and the assumptions become applicable with 
A = a" 1 and B = a~ 2 . The correlation functions are the same both in the cases of the 
usual and renormalized distributions. If we use the renormalized distributions and other 
assumptions, we can derive Eq. ( |B12| ), but the derivation is more complicated. 
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FIGURES 





<N> 

FIG. 2. Chaoticity and weight factor as functions of the mean number of coherent sources in 
the model of multi-coherent sources. Solid and dashed lines stand for the chaoticity and weight 
factor, respectively. 
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FIG. 4. Weight factor as a function of a. The model and lines are the same as in Fig. ||. 
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FIG. 5. Weight factors as functions of chaoticities. Solid and dashed lines stand for the par- 
tially coherent source and the multi-coherent sources, respectively. The plot with error-bars is the 
experimental data from NA44. 




FIG. 6. Weight factors as functions of chaoticities in the model of multi-coherent sources and 
one chaotic source, varying e from to 1. The lines from down to up stand for a = 0.1, 0.2, 0.4, 
0.6, 0.8, 1.0, 1.5, 2.0, respectively. The plot with error-bars is the experimental data from NA44. 
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